Chaos-assisted emission from asymmetric resonant cavity microlasers 
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We study emission from quasi-one-dimensional modes of an asymmetric resonant cavity that are 
associated with a stable periodic ray orbit confined inside the cavity by total internal reflection. 
It is numerically demonstrated that such modes exhibit directional emission, which is explained 
by chaos- assisted emission induced by dynamical tunneling. Fabricating semiconductor microlasers 
with the asymmetric resonant cavity, we experimentally demonstrate the selective excitation of the 
quasi-one-dimensional modes by employing the device structure to preferentially inject currents to 
these modes and observe directional emission in good accordance with the theoretical prediction 
based on chaos-assisted emission. 

PACS numbers: 42.55.Sa, 05.45.Mt, 42.60.Da, 42.55.Px 



U 



> 



o 



X 



I. INTRODUCTION 

Besides growing interest from the fundamental physics 
of light and practical applications in photonics , two- 
dimensional optical microcavities have attracted much 
attention from the field of quantum/ wave chaos, since 
they can be regarded as physical realization of dynami- 
cal billiards, a well-studied model in the field (for reviews 
see for example Refs. 0-01 )• Concepts of the quan- 
tum/wave chaos theory have been successfully applied to 
two-dimensional microcavities to explain how ray chaos 
influences the characteristics of resonant modes such as 
emission directionality and lifetimes. 

One of such applied concepts is dynamical tunneling 
Q, which refers to tunneling between classical invariant 
components disconnected from each other under classical 
dynamics (see Ref . @ for a review) . In generic Hamilto- 
nian systems, i.e., those with mixed phase space, where 
regular and chaotic dynamics coexist, dynamical tunnel- 
ing manifests itself in the fundamental processes such as 
reg ular-to-chaotic and chaotic-to-regular tunneling fl^ 
|l3|. One of the dramatic effects of dynamical tunneling 
is observed in the enhancement of level splitting of dou- 
blets associated with symmetry-related regular compo- 
nents separated by a global chaotic sea. The phenomenon 
is explained by the mechanism known as chaos-assisted 
tunneling 13" 13) the regular components are coupled 
not only by direct regular-to-regular tunneling but also 
by a multi-step process involving regular-to-chaotic and 
chaotic-to-regular tunneling, and this process leads to a 
significant enhancement of the coupling. The dynamical 
tunneling and chaos-assisted tunneling have been experi- 
mentally observed in microwave billiards [20,., 21.] and cold 
atom systems [l^ [l^ . 

In the context of optical microcavities, the effect of 
dynamical tunneling was already pointed out in the pio- 
neering theoretical study of asymmetric resonant cavities 
[23 |. where the significant shortening of resonance life- 



times is attributed to dynamical tunneling. This effect 
was quantitatively studied by effective Hamiltonian mod- 
els |25l - [27| , where a formula is derived describing how res- 
onance lifetimes depend on the factors such as the area 
of regular components and wavenumber. Strong influ- 
ence of dynamical tunneling on the emission pattern of a 
ray-chaotic cavity was also pointed out in Refs. (26l[27|. 
Resonance lifetimes were connected with a direct regular- 
to-chaotic tunneling rate in Ref. f28j . The enhancement 
and fluctuation of frequency splitting of regular mode 
doublets is explained by the chaos- assisted tunneling in 
Ref. [29]. 

Recently, the manifestation of dynamical tunneling 
was experimentally observed in the emission patterns 
of semiconductor microlasers [s^, HH and the pump- 
coupling efficiency of dye-doped liquid jet microlasers 
(3^ . In the former study, an asymmetric resonant cav- 
ity was designed to have quasi-one-dimensional modes 
strongly localized along a stable periodic ray orbit that 
is confined by total internal reflection (TIR). In the ex- 
periment, such modes were excited by selective pumping, 
and highly directional emission was observed. The pur- 
pose of this paper is to present a comprehensive study 
on this microlaser. The emission directionality of this 
microlaser was interpreted by the following two steps: 
(i) dynamical tunneling from the TIR-confined ray orbit 
to a neighboring chaotic ray orbit, and (ii) ray-dynamical 
propagation following the chaotic dynamics leading to re- 
fractive emission. In this paper, we shall call this process 
chaos-assisted emission in order to emphasize that the 
process is regular-to-chaotic-to-continuum. 

This paper is organized as follows: In Sec. |TT1 we 
introduce the asymmetric resonant cavity and explain 
its ray dynamics, where particular attention is paid to 
the rectangle-shaped stable ray orbit confined by TIR. 
Section IIIII presents the numerical study of the resonant 
modes associated with this orbit (rectangle-orbit modes, 
in short), showing that these quasi-one-dimensional 
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modes exhibit chaos-assisted directional emission. Sec- 
tion IIVI explains the device structure of fabricated semi- 
conductor microlasers, especially focusing on the struc- 
ture to achieve selective excitation of the rectangle-orbit 
modes. Section |V] provides measured spectral data that 
evidence the selective excitation. Section IVII presents 
measured emission patterns which are attributed to the 
chaos-assisted emission. Finally, our conclusion is given 
in Sect. EIIl 



II. ASYMMETRIC RESONANT CAVITY 

In this paper, we consider the asymmetric resonant 
cavity whose shape is defined in the polar coordinates 
{r,(f)) as 

r — R[(f)) = ro(l -I- acos2(/) -f & cos 40 -I- ccos6(/)), (1) 

where is the size parameter and a, 6, and c are deforma- 
tion parameters fixed as a = 0.1, h — 0.01, and c — 0.012. 
The geometry of the cavity is depicted in Fig. [T] (a) . 

An important feature of this cavity is the existence of 
the 4-bounce stable ray orbit drawing the "rectangle" as 
shown in Fig. [I](b), which is confined by TIR when the 
refractive index n > ^/2. The peculiarity of the rectangu- 
lar orbit in relation with global ray dynamics is revealed 
by plotting the Poincare surface of section (PSOS) of 
the ray dynamics. The PSOS describes the successive 
bounces of a ray trajectory usually using the Birkhoff 
coordinates (s,sin0), where s represents the arclength 
along the cavity and the incident angle as illustrated 
in Fig. [T](a). The PSOS is convenient for describing the 
openness of the dielectric cavity, because the condition 
for TIR is given by |sin0| > 1/n. Throughout this pa- 
per, we assume n = 3.3. This value corresponds to the 
effective refractive index calculated for the structure of 
actually fabricated devices to be explained in Sect. IIVI 

Because the PSOS is symmetric under sin 61 i~?> — sin6', 
we show only the upper half of the PSOS in Fig. [1] (c), 
where one finds that ray dynamics is mostly chaotic, and 
the four islands of stability corresponding to the rect- 
angular orbit are the only dominant islands located in 
the TIR-region, i.e., sm9 > 1/3.3. In this paper, we fo- 
cus on the resonant modes that are strongly localized on 
the TIR-confined dominant islands. The ray-dynamical 
property that these islands are isolated and only dom- 
inant islands in the global chaotic sea will later enable 
us to interpret the emission patterns of those resonant 
modes by a simple scenario based on chaos- assisted emis- 
sion. 



III. RECTANGLE-ORBIT MODES 

Resonant modes of the cavity are given as the solu- 
tions of the Helmholtz equation derived from Maxwell's 
equations [1, H, H^l . They are quasibound states with fi- 
nite lifetimes characterized by complex wavenumbers k. 
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FIG. 1: (color online) (a) The geometry of the asymmetric 
resonant cavity with the definitions of the polar coordinates 
(r, (/!>) and of the Birkhoff coordinates (s,sin0). (b) The rect- 
angular ray orbit, (c) The upper half of the Poincare surface 
of section (PSOS) for the ray dynamics. The coordinate s is 
normalized by the perimeter of the cavity. The critical line 
for total internal reflection, sinO — 1/3.3, is plotted by a dot- 
ted line. The four flUed squares (■) represent the rectangular 
orbit shown in (b). 



In this paper, we are interested in quasi-one-dimensional 
modes strongly localized along the stable rectangular or- 
bit. The existence of such modes can be proven in the 
short-wavelength limit by the so-called Gaussian optics 
[2^. This type of modes has been previously studied, 
for instance, for the so-called 'bow-tie' laser (ssj, but 
an important difference here is that the corresponding 
ray orbit is confined by TIR. Thus, the ray optics based 
on Fresnel's law as well as the Gaussian optics predict 
that rectangle-orbit modes are perfectly confined, which 
is not true for accurate wave description. In order to 
correctly capture the emission patterns of the rectangle- 
orbit modes, below we resort to numerics. 

An accurate way for numerically calculating the modes 
is provided, for instance, by the boundary element 
method [s^l • In Fig. [2] (a) we present a result obtained 
by this method, showing the field intensity distribution 
of a resonant mode strongly localized along the rectangle 
orbit, where one can see a regular modal pattern along 
the orbit. Its wavenumber is A: = 49.94 - i 0.00012. Con- 
cerning the resonant mode calculation, we note that the 
average radius of the cavity is fixed as tq = 1. We em- 
ploy the boundary condition for the transverse-electric 
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polarization, to be consistent with experimental data pre- 
sented later. 

In order to observe exponentially small emission from 
the rectangle-orbit mode, we plot the field intensity dis- 
tribution in log scale in Fig. [2] (b), where one finds that 
dominant emission occurs on both sides of the cavity to- 
wards the directions (/) = ±90°. These emissions are un- 
derstood as the wave tunneling effect intrinsic to a ray- 
chaotic cavity. This emission mechanism is better ex- 
plained by the PSOS than the real space. Thus, we con- 
sider the PSOS-representation of an eigenfunction called 
the Husimi distribution [H, 111] . Figure [3] (a) shows the 
Husimi distribution of the rectangle-orbit mode shown in 
Fig. [2] (a), where it is plotted in log scale to see expo- 
nentially small intensities in the area of the chaotic sea. 
In Fig. [21(a), the intensity increases as the color changes 
from black to white. Refiecting the fact that the mode is 
localized along the rectangular orbit, we find high inten- 
sities on the corresponding islands. In addition, we find 
some intensities spread over the chaotic sea, even reach- 
ing down to the critical lines for TIR, sin 6* = ±1/3.3. 
The four high-intensity spots on the critical line labeled 
a, b, c, and d correspond to the four dominant emissions 
observed in Fig. [2] (b) . Together with the field intensity 
distribution in the real space (Fig. [2] (b)), the Husimi 
distribution reveals that emission occurs at two close but 
different boundary points (one towards 4> = 90° and the 
other towards = —90°) on each side of the cavity. 

The emission pattern of the rectangle-orbit mode 
can be qualitatively explained in the following manner 
[26l [27| . Although intensities are highly localized on 
the islands corresponding to the rectangular orbit, some 
intensities are leaked to their neighboring chaotic or- 
bits via dynamical tunneling. Once this happens, the 
leaked intensities are transported to the leaky region (i.e., 
I sin^l < 1/3.3) by the chaotic dynamics, resulting in re- 
fractive emission. Hence, it is the ray dynamics around 
the critical lines for TIR that determines the emission 
pattern. In fact, such ray dynamics is governed by the 
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FIG. 2: (color online) (a) Field intensity distribution of the 
rectangle-orbit mode with the wavenumber k = 49.94 — 
j 0.00012. (b) Log-scale plot of the rectangle-orbit mode. The 
intensity increases as the color changes from black to white. 



unstable manifolds emanating from a short-periodic un- 
stable periodic point located near the critical line 
l4ll |. In Fig. |3] (b), we show how a ray ensemble initially 
prepared around the leftmost island is mapped to after 
4 bounces. The initially circular area is stretched and 
folded after 4 bounces, partially penetrating into the re- 
gion below the critical line as magnified in the inset of 
Fig. [21(b). As iUustrated in Fig. [21(a), the high-intensity 
spot labeled c coincides with the position where the un- 
stable manifold crosses the critical line. The positions of 




FIG. 3: (color online) (a) Log-scale plot of the Husimi distri- 
bution for the rectangle-orbit mode with k = 49.94— i 0.00012, 
where the islands corresponding to the rectangular orbit 
and the unstable manifolds of the unstable period-3 point 
(0.5,0.389...) are superimposed. The intensity increases as 
the color changes from black to white. The critical lines for 
total internal reflection, sinO = ±1/3.3, are plotted by dotted 
lines, (b) Ray-dynamical evolution of an ensemble initially 
prepared over the leftmost island and its neighborhood in the 
upper half PSOS. The inset shows a magnification around the 
area where the unstable manifold overlaps with the critical 
line. 
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the other high-intensity spots (a, b and d) are explained 
in the same manner. 

In the remaining part of this paper, we study whether 
the theoreticahy predicted chaos-assisted emission from 
the rectangle-orbit modes can be experimentally ob- 
served. Therefore, it is crucial in experiments to excite 
only the rectangle-orbit modes, suppressing the excita- 
tion of the other modes. In order to achieve this selective 
excitation, we employ the selective pumping method that 
has been previously examined for quasi-stadium-shaped 
microcavity lasers, where the selective excitation of two 
different types of modes, called the axis modes and ring 
modes, has been successfully demonstrated (43 - l48| . 



test the sensitivity of the selective excitation on a contact 
window pattern, we examine three different patterns of 
the contact window, i.e., a broad, an intermediate, and 
a narrow one, respectively shown in Figs. [5] (a), [5] (b), 
and[5](c). We show in Fig. [5](e) how the narrow contact 
window overlaps with the rectangle-orbit mode. 

On the rim of the cavity, we make the "margin area" 
where Si02 is left. This is to avoid excitation of 
whispering-gallery-type modes that are localized along 
the rim of the cavity. The precise shape of the margin 
area is given by 0.75i?((/)) < r < R{(t>)- The p-electrode 
metal is deposited over the contact window and part of 



IV. DEVICE STRUCTURE 

In this section, we explain the device structure to 
achieve the selective excitation (see Ref. [i^ for the 
details of the fabrication process). Figure |4] shows a 
schematic cross-sectional view of the device. The de- 
vice is fabricated from MOCVD-grown material with a 
1.5-/Ltm n-Alo.5Gao.5As lower cladding layer, a 0.2-/im 
n-Al2;Gai_a;As {x = 0.5 — 0.2) graded region, a 10-nni 
GaAs quantum well layer, a 0.2-/im p-A\xGa.i-x^s{x = 
0.2 — 0.5) graded region, a 1.5-/im p-Alo.5Gao.5As up- 
per cladding layer, a 0.2-//m p-GaAs contact layer, and a 
400-nm Si02 layer. The average radius ro of the device 
is fixed as 50 /xm. 

The key structure for the selective excitation is the 
"contact window" etched to the Si02 layer. The Si02 
layer works as the insulation between the semiconductor 
layers and the p-clcctrode metal, enabling a selective cur- 
rent injection through the contact window. The top-view 
shape of the p-GaAs layer is given by adding l-/.tm mar- 
gins to the regions of the contact window which is formed 
by a dry etching process before depositing the Si02 layer. 

In order to excite only the rectangle-orbit modes, we 
form the contact window along the rectangular orbit as 
illustrated in Figs. [S](a)-(c). It is not known a priori what 
design of the contact window is most optimal for stable 
and efficient excitation of the rectangle-orbit modes. To 
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FIG. 4: (color online) Schematic cross-sectional view of the 
device. 
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FIG. 5: (color online) (a)-(c): Top-view geometry of the Si02 
layer. The white regions inside the cavity define the contact 
window where Si02 is removed and through which currents 
are injected. The contact window is broad in (a), intermediate 
in (b), and narrow in (c). We make the "margin area" on the 
rim, whose shape is defined by O.75-R(0) < r < R(<j)). Above 
the Si02 layer, the p-electrode metal is deposited as in (d), 
whose shape is defined by r < 0.8i?(<^). (e) Field intensity 
distribution of the rectangle-orbit mode with k = 49.94 — 
iO. 00012 superimposed with the narrow contact window, (f) 
Scanning electron microscope image of a fabricated device 
with the broad contact window. The average radius is ro = 50 
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the surrounding Si02 layer as shown in Fig. [5](d) and its 
shape is defined by r < 0.8R{4>). In Fig. [5] (f), we show 
the scanning electron microscope image of a fabricated 
device with the broad contact window. 



V. SPECTRAL CHARACTERISTICS 

As, for instance, the Gaussian optic theory directly 
demonstrates [ll] , rectangle-orbit modes are regarded as 
the quantization along the rectangular orbit. Their lon- 
gitudinal mode spacing is given by 
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FIG. 6: (color online) Lasing spectra of the device with the 
narrow contact window for different pumping currents: (a) 50 
mA; (b) 80 mA; (c) 100 mA. Equidistantly spaced peaks are 
observed in (c), where the average mode spacing is 0.67 nm, 
in good agreement with the longitudinal mode spacing of the 
rectangle-orbit modes, 0.63 nm. 




where A is the wavelength, n the effective refractive in- 
dex, L the path length of the rectangle orbit, and dn/dX 
the dispersion. In this section, we present measured 
spectral data showing equidistantly spaced peaks corre- 
sponding to the above mode spacing, which constitutes 
definitive evidence that rectangle-orbit modes only are 
excited. The devices are tested in pulsed mode opera- 
tion at 25°C with pulse width 500 ns and repetition rate 
1 kHz. In Fig. [6l we show lasing spectra of the device 
with the narrow contact window for the pumping cur- 
rents 50 mA, 80 mA, and 100 mA. For 50 mA, which 
is just above the lasing threshold current, one can con- 
firm single-mode lasing. For 100 mA, the spectrum shows 
equidistantly spaced peaks with the averaged mode spac- 
ing AA = 0.67 nm. Considering the uncertainty of the 
refractive index and the resolution of the measurement, 
one can conclude that this value corresponds to the lon- 
gitudinal mode spacing of the rectangle-orbit modes esti- 
mated as 0.63 nm. This theoretical estimate is calculated 
from Eq. © with A = 860 nm, n = 3.3, L = 283.04 fim, 
and dn/dX = —1.0 fj,m~^. Here we use the value of the 
dispersion dn/dX from Ref. poj . 

Because currents are not injected around the rim of 
the cavity, the lasing of whispering-gallery-type modes is 
highly unlikely. This can be confirmed by estimating the 
mode spacing corresponding to whispering-gallery-type 
modes. Using the above formula with L — 317.75 /xm 
(i.e., perimeter of the cavity), one can estimate their 
mode spacing as 0.56 nm, which is considerably smaller 
than the measured mode spacing. 

The above theoretical estimates convince us that the 
equidistantly spaced peaks in Fig. [6](c) are the rectangle- 
orbit modes, not the whispering-gallery-type modes. In 
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FIG. 7: (color online) Lasing spectrum of the device with 
the narrow contact window for the pumping current 500 mA. 
The inset shows the magnification of a high-intensity range, 
where equidistantly spaced peaks are observed, whose average 
mode spacing is 0.64 nm, in excellent agreement with the 
longitudinal mode spacing of the rectangle-orbit modes, 0.63 
nm. 
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addition, using Fig. |6] (c) as a reference, we can also 
conclude that peaks observed in Figs. H] (a) and (b) cor- 
respond to the rectangle-orbit modes. 

We confirmed selective excitation of the rectangle-orbit 
modes for a much higher pumping current. In Fig. [71 we 
show the lasing spectra for the pumping current 500 mA, 
where one finds again equidistantly spaced peaks with 
the average mode spacing 0.64 nm, which is in excellent 
agreement with the theoretical estimate 0.63 nm. The 
spectral peaks in Fig. [7] are red-shifted from those in 
Fig. [HI which is considered to be due to a heat effect 
caused by the high pumping current. 

Concerning the lasing of higher transverse excited 
modes, we cannot find its clear evidence for lower pump- 
ings. However, for the pumping 500 mA, we can identify 
subpeaks at, for instance, 861.7 nm and 862.4 nm, which 
are considered to be higher transverse excitations. 

We also confirmed that the selective excitation is 
achieved with the intermediate and broad contact win- 
dows. Figures m (a) and|8](b) show spectra for the device 
with the intermediate contact window at 140 mA and for 
the device with the broad contact window at 80 mA, re- 
spectively. In both cases, one can confirm equidistantly 
spaced peaks. The average mode spacing is 0.67 nm in 

3000 I ■ ■ ^ ■ ■ ' 

. (a) Intermediate, j 
140 mA 
2000 ■ I 

i Hi 




' ' ' ' ' ' ' ' 

850 852 854 856 858 860 862 864 
Wavelength (nm) 



FIG. 8: (color online) Lasing spectra for (a) the device with 
the intermediate contact window at the pumping current 140 
mA and for (b) the device with the broad contact window 
at 80 mA. Equidistantly spaced peaks are observed in both 
cases. The average mode spacing is 0.67 nm in (a), while 
0.68 nm in (b), both in good agreement with the longitudinal 
mode spacing of the rectangle-orbit modes, 0.63 nm. 



Fig. [8] (a) while 0.68 nm in Fig. |8] (b). These values are 
again in good agreement with the theoretical estimate 
0.63 nm. As for the lasing threshold currents, we could 
not find significant dependence on the contact window 
patterns, i.e., the threshold currents are around 50 mA 
irrespective of the contact window pattern. 



VI. EMISSION CHARACTERISTICS 

For microcavities with mostly or fully chaotic ray dy- 
namics, it has been reported that experimental far-field 
patterns are well explained by ray-dynamical calculations 
t4Q . i43l . l50l - [53| . Generally, some intrinsic deviation is ob- 
served between the far-field pattern of a single resonant 
mode and that of ray-dynamical calculation, especially 
in the short-wavelength regime [51 54]. Nonetheless, it 
was demonstrated that the averaged far-field pattern of 
many resonant modes is very well reproduced by ray- 
dynamical calculation. Because multi-mode lasing works 
as a mechanism to realize the "average" of many modes, 
one can expect that the correspondence between experi- 
mental data and ray-dynamical calculation becomes the 
better, the more the number of lasing modes increases. 
In fact, the validity of this expectation was confirmed for 
a stadium-shaped microlaser [s^l- In this section, we re- 
port measured emission patterns of our devices where the 
rectangle-orbit modes are selectively excited and show to 
which extent they can be reproduced by ray-dynamical 
calculation. 

We show far-field patterns of the device with the nar- 
row contact window for the pumping currents 100 mA 
and 500 mA in Figs. H] (a) and[H](b), respectively. We 
note that these patterns are normalized after subtract- 
ing uniform background due to spontaneous emission. In 
good accordance with the theoretical prediction based on 
the chaos-assisted emission, we observe strong emissions 
towards ±90°. 

In Figs. M (a) and [9] (b), we superpose the theoretical 
far-field pattern obtained by ray-dynamical calculation. 
We find that the experimental far-field pattern for 500 
mA is well reproduced by ray-dynamical calculation even 
for the substructures of the main peaks. As shown in Ref. 
[sO], wave calculation of the far- field pattern for a single 
mode does not show such substructures (see Fig. 5 of 
Rcf. ^3(7|). Therefore we conclude that the substructures 
of the main peaks in the measured far-field pattern are 
due to multimode lasing, because the far-field pattern for 
ray-dynamical calculation can be considered as that from 
an average of many modes as explained above. 

For 100 mA, the experimental far-field pattern is asym- 
metric, showing some deviations from the ray-dynamical 
calculation. An intrinsic mechanism leading to the asym- 
metry is the locking of modes belonging to different par- 
ities [1^, [5^ . The locking is caused by nonlinear modal 
interaction due to the lasing medium, creating a lasing 
mode whose emission pattern violates the symmetries of 
the cavity. When the number of such a lasing mode 
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FIG. 9: (color online) Far-field emission patterns of tiie device 
with the narrow contact window for (a) 100 mA and (b) 500 
mA. The definition of the angle </!> is given in Fig. [l](a). The 
solid curves represent experimental data, while the dotted 
curves represent ray-dynamical data. 



is large, the asymmetry is averaged out, resulting in a 
symmetric emission. However, the asymmetry becomes 
prominent when the number of such a mode is relatively 
small ^55i]. For the pumping 500 mA, we can identify 
many modes with significant intensities in the lasing spec- 
trum shown in Fig. 7, while for the pumping 100 mA we 
find only a couple of dominant lasing modes in the las- 
ing spectrum shown in Fig. 6 (c). We expect that for 
100 mA the relatively small number of the lasing modes 



FIG. 10: (a) CCD photo of the lasing device with the broad 
contact window for 300 mA. The photo is taken from a slightly 
oblique angle, causing the difference between the right hand 
edge and left hand edge. The outline of a needle for current 
injection is indicated by the dotted line, (b) The field inten- 
sity distribution of the rectangle-orbit mode shaded by the 
electrode metal area. 



makes the emission asymmetry significant. 

Lastly, we show a CCD photo of the lasing device with 
the broad contact window at 300 mA in Fig. [10] (a). 
Although the lasing emission is in-plane, some light scat- 
tered inside the cavity is observed from the top in the 
margin area. For comparison, we show in Fig. [10] (b) the 
field intensity distribution of the rectangle-orbit mode 
shaded by the electrode metal area. We can see that 
bright regions in Fig. [10] (a) are coming mainly from 
the four corners of the rectangular orbit. Moreover, we 
can observe two bright spots around the right hand edge 
of the cavity. These are considered to be the light scat- 
tered at the cavity boundary, thus indicating the emission 
points. These emission points are in good accordance 
with the theoretical prediction drawn from the Husimi 
distribution shown in Fig. [3] (a) as well as the near-field 
pattern shown in Fig. [2] (b). 



VII. CONCLUSION 

We have focused on the quasi-one-dimensional modes 
of the asymmetric resonant cavity, which are strongly 
localized along the stable rectangular ray orbit that is 
confined by TIR. Direct numerical calculation of such a 
resonant mode revealed that directional emission occurs 
from cavity boundary points apart from the rectangu- 
lar ray orbit. This at first sight counterintuitive result 
was explained by the chaos-assisted emission mechanism 
consisting of regular-to-chaotic tunneling and subsequent 
chaotic transport leading to refractive emission. 

We fabricated the cavity with the structure to selec- 
tively excite the rectangle-orbit modes. The success of 
the selective excitation was confirmed by the observa- 
tion of the equidistantly spaced peaks in the lasing spec- 
tra well corresponding the longitudinal mode spacing of 
the rectangle-orbit modes. In the emission patterns, we 
observed directional emission attributed to the chaos- 
assisted emission. 

Our results demonstrated that the concepts of quan- 
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turn/ wave chaos and nonlinear dynamics such as dynam- 
ical tunneling and unstable manifolds are indispensable 
for precisely understanding phenomena occurring in the 
micro devices. As a way to extract light from optical cav- 
ities, the chaos-assisted directional emission adds a new 
idea to the conventional theories [HH , and it would be of 
interest to apply this idea to design practical devices. 
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